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We consider a system of bosonic atoms in an axially symmetric harmonic trap augmented with a two dimen- 
sional repulsive Gaussian optical potential. We find an expression for the grand free energy of the system for 
configurations ranging from the harmonic trap to the toroidal regime. For large tori we identify an accessible 
regime where the ideal gas thermodynamics of the system are found to be independent of toroidal radius. This 
property is a consequence of an invariant extensive volume of the system that we identify analytically in the 
regime where the toroidal potential is radially harmonic. In considering corrections to the scale invariant transi- 
tion temperature, we find that the first order interaction shift is the dominant effect in the thermodynamic limit, 
and is also scale invariant. We also consider adiabatic loading from the harmonic to toroidal trap configuration, 
which we show to have only a small effect on the condensate fraction of the ideal gas, indicating that loading 
into the scale invariant regime may be experimentally practical. 
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I. INTRODUCTION 

While Bose-Einstein condensation of dilute gases is now 
routinely observed, the degree of control and interrogation 
now affords more detailed studies of the dynamics of the con- 
densation process, such as the symmetry breakingassociated 
with the Kibble-Zurek mechanism (KZM) fl, Bll 0- Re- 
cent experiments loading toroidal traps with Bose-degenerate 
gases Jj, @] have shown the need for a basic theoretical un- 
derstanding of the properties of Bose-Einstein condensation 
in toroidal traps, which may be important for future studies of 
BEC formation and tests of KZM. 

Following the development of storage rings for neutral 
atoms 01, toroidal traps for BECs were proposed usingvari- 
ous combinations of magnetic and optical techniques IH @] . 
Sagnac interferometry is an important application for large 
toroidal traps, having the feature of resolution proportional 
to the area of the interferometer fioll . While traps of order 
~ 3 mm diameter have been created fill [l2tl . BEC loading 
for such large traps was limited to launching the BEC into 
the toroid which then acts as a dispersive waveguide. Smaller 
toroids, which can be more easily loaded, have recently been 
produced JH 0, [H, [3 opening the way for studies of super- 
fluidity and persistent currents in a nontrivial trapping topol- 
ogy- 
Theoretical efforts have focussed on BECs far below T c , 
Many features of toroidal geo metry have been studied, in- 
cluding topological phases lfl5l|. th e stability of macroscopic 
persistent currents flril ITtL Il8lll9l I20I l2lll , excitation spec- 
tra l22tl . ato mic p hase interference devices l23ll . vortex- vortex 
interactions 12411 . generation of excitations via stirring l25ll . 
dynamics of sonic horizons l26ll . parametric amplification of 
phonons IT27ll . rotational current generation l28tl . the inter- 
play of interactions and rotation IT29ll . giant vortices l30ll . and 
vortex signatures l3lll . Ideal gas theory has recently been 
used ll32ll to study the rapidly rotating Bose gas in a quartically 
stabilized harmonic trap realized at ENS 13311 . The BEC tran- 
sition temperature in non-power law traps is thus becoming 
more relevant, and a recent study of optical lattices [34] gives 
further indication that analytical expressions can be found for 



increasingly rich potentials. However, the ideal gas thermo- 
dynamics of three dimensional toroidal potentials-crucial for 
understanding the dynamics of Bose-Einstein condensation- 
have not been addressed. 

The effect of trapping geometry on the BEC transition was 
emphasized by Bagnato et al. Il35ll . It was observed that, for 
power law traps, increasing the confinement of the system has 
the effect of increasing the peak phase space density and thus 
raising the temperature of the BEC transition. An understand- 
ing of how phase space density depends on the toroid size is 
crucial for making large ultra-cold persistent currents. There 
is also the role of topology to consider. In particular, in a 
toroidal trap the angular spatial coordinate becomes unavail- 
able for thermalization, suggesting a potentially interesting in- 
terplay between topology and system size. 

In this work we study the properties of a Bose gas trapped 
by a specific toroidal potential. The potential is created from a 
harmonic magnetic potential combined with a repulsive (blue 
detuned) optical potential with a Gaussian spatial profile Ja); 
we refer to this potential as harmonic-Gaussian and show that 
it has uniquely interesting properties which advantage it for 
creating large toroidally trapped BECs. Using the semiclassi- 
cal approach to the thermodynamics of the ideal gas we find 
an exact expression for the free energy. 

Examining the properties of the system for increasing 
toroidal radius shows the existence of an analytically tractable 
regime of scale invariance with respect to the toroid radius. In 
this regime the toroidal trap is well approximated as radially 
harmonic. We use the theory of Romero-Rochin l36ll to iden- 
tify the invariant generalized extensive volume of the system 
which governs scale invariance. We further generalize this re- 
sult to show that the system enters a scale invariant regime 
even when this approximation is not valid. Focusing on the 
preservation of quantum degeneracy, we then treat the BEC 
transition temperature in detail and consider finite size and 
mean field corrections to the scale invariant result. Finally, we 
discuss possible means to reach the scale invariant regime. 
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FIG. 1 : The harmomic-Gaussian trapping potential is shown at z = 
(solid lines) for an optical potential of width cr = 70pm and (a) 
V /k B = 124 nK, (b) V /k B = 211 nK, (c) V /k B = 557 nK, and 
(d) Vo/k B = 2784 nK. In (a) the dashed-dotted curve is the bare 
harmonic potential (Vo = 0) and the dotted line is V a = mtofcri/2, 
which in this case is also the critical case for bifurcation of the trap 
minimum: Vq = Va-. The dashed lines in (b)-(d) give the harmonic 
approximation to the radial potential. The harmonic trap frequencies 
are (to z ,(o r ) = 2^(15.3,7.8) Hz 



A. Geometry of the harmonic-Gaussian potential 

We consider a Bose gas confined in the trapping potential 
V(x) = ^ (co 2 r r 2 + to 2 z 2 ) + Vop(r), (1) 



where axial and radial trapping frequencies are co z , at r , and the 
non-harmonic potential is given by 



V 0P (r) = Voexpr-rVoX 



(2) 



where r is the distance from the z-axis. This potential can 
be created using a magnetic trap combined with a detuned 
laser field propagating along the z-axis which forms an op- 
tical dipole potential 13711 . To determine the geometry of the 
trap it is convenient to define the energy 



Va- = ^mco 2 r crl, 



(3) 



which is the potential energy of an atom at r = <x in the 
harmonic potential; this energy will play a central role in de- 
termining the thermodynamics of the system. The minimum 
of the combined potential is located at z — and r - r,„, where 



[o for v <v a , 

UcrllniVo/Vcr) for V > V v , 



(4) 



from which we find that V m = min[V(x)] = V(r = r m ,z = 0) 
is given by 



JVo for Vb<V«-. (5) 

'WM+MVo/V,,)) for V () >V LT . 



There are three distinct regimes parametrized by the ratio 
Vo/V lT which are 

1. Dimple trap {Vq < 0) The Gaussian forms a dimple in 
the center of the harmonic trap. 

2. Flat trap (0 < Vo < Vo-) The trap is flattened but not 
toroidal, as shown in Fig.[T](a). 

3. Toroidal trap {V a < Vq) The trap becomes toroidal, as 
shown in Fig.[T](b)-(d). 



H. IDEAL BOSE GAS IN A HARMONIC-GAUSSIAN TRAP 

In this section we develop the general grand-canonical the- 
ory of the ideal Bose gas in the harmonic-Gaussian trap. 

A. Grand canonical free energy 

In general, the grand potential function for the system is 

r = -k B T In Z, (6) 

where Z is the grand partition function. For the Bose gas 
distributed over levels with excitation energy e, this becomes 



;F = fc B r^ln(l -e^- e,) ), 



(7) 



f3 = 1/kgT. For the situation of interest, the chemical potential 
is assumed to approach the ground state energy of the trap 
eo, leading to macroscopic occupation of the ground state A^. 
In this regime, using the semi-classical approximation for the 
excited states, T can be written as 



/ m \3/2 



\27Tfc 2 



rz 



(k/3) 5 l 2 



where 



Girt 



-s 



d x exp [-o'V(x)]. 



om. (8) 



(9) 



All thermodynamic properties are determined by @(a) and its 
derivatives. The BEC transition temperature is found from the 
total atom number at the point where the chemical potential 
reaches the ground state energy 



dp 



(10) 
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with jVo = 0, giving the transition temperature T c as the solu- 
tion of 



m, 



NA 



(ID 



where Ajb = ^2nh 2 IrnksT is the thermal de Broglie wave- 
length. 

B. Thermodynamics of the harmonic-Gaussian trap 

For the general form of the harmonic-Gaussian trap ([]]) we 
can evaluate (O to find 



2 I n yiaVo-, aV ) 



(aV ) a 



where 



y(a, x) 



r e-'r 

Jo 



dt. 



(12) 



(13) 



is the incomplete Gamma function. 

In the semiclassical approximation it is consistent to take 
the ground state energy as the trap minimum eo — > V m , and 
the free energy can now be written as 



r = N (.V m -fi) 



/3 4 tt 3 cd-oj, 

where we define the generalized ^"-function 



(14) 



f v (z, a,b) = a 2^ -j—[^(ka)y*(ka, kb), (15) 

and j* (a, x) = x~ a y(a, x)IT{a) is a sin gle v alued analytic func- 
tion of a and x with no singularities 13811 . This function has 
the following asymptotics: 

Km&CeW.^jffVb) = (16) 
lim ° ^f^°> = (17) 



lim f v («*',j8V'<r,^Vo)^ = W^™'), (18) 

where ( y (z) = 2^=o z */^ V ' s trie polylogarithm function. In 
the regime of Bose-Einstein condensation, where fj. - V,„ — > 
this further reduces to the ordinary Riemann-zeta function 

lim^o-TvCl./^ySVo) = = ?(v). 

The number of atoms in the system is given by 



N ~ ~ N ° + oW 2 ' 

dfi p 3 n s a}$tt) z 



(19) 



from which the transition temperature T c for N atoms is found 
as the solution of 



N = 



P 3 c ti 3 a>}(o z 



(20) 



where f} c = l/kBT c . From equations ( fT9l ) and ( f20b we then 
find the equation of state as 



N \T r 



h(ee v »,pv<r,pv ) 



(21) 



which must be evaluated numerically and describes all the 
three regimes introduced above: dimple trap, flat trap, and 
the toroidal trap. 

We remark that all ideal gas properties of the system for ar- 
bitrary ctq and Vo are not easily obtained from the current for- 
mulation. The difficulty arises when considering derivatives 
of ([TBI l with respect to a and b, as is required for obtaining 
the entropy and heat capacity. The function (TT3T > is not closed 
with respect to differentiation, rather leading to a hierarchy of 
transcendental functions with each successive derivative oper- 
ation lf38h . This problem can be solved by introducing a gener- 
alization of ([TBI l. and a treatment that includes the double well, 
toroidal and ellipsoidal cases will be provided elsewhere. For 
the remainder of this paper we will restrict our attention to the 
< Vo case and focus primarily on the toroidal regime. 



III. HARMONIC SCALE INVARIANCE 

So far we have seen that known results are obtained in the 
appropriate limits. Our new result is found by using Eq. ( TP71 ), 
where we find for the toroidal trap 



r = N (V m -fi)- 



V27TCr £ 7 / 2 ( 



AdB /3 3 h 2 u) r u) z 



(22) 



What is immediately apparent here is that in this regime all 
of the ideal gas thermodynamical properties become indepen- 
dent of Vo. In particular, the thermodynamics are independent 
of the toroidal size, and the system enters a scale invariant 
regime defined by 



k B T « V a « V . 



(23) 



In this regime both of the energy scales Vo and V a drop out 
of the problem. Specifically, for fixed harmonic frequencies 
(w r , (o z ), and a beam of fixed width (cro), the laser power 
can be increased to generate a toroidal potential with larger 
perimeter. Nevertheless, all ideal gas properties of the sys- 
tem are invariant under this dilation. The only remnant of the 
Gaussian beam enters through the length scale <t appearing 
in Eq. d22l . 

The invariance condition d23l can also be written in terms 
of appropriate length scales of the trap. We introduce a tem- 
perature associated with the optical potential height Vo = 
£bJo, and de Broglie wavelength Aq = ^2nh 2 IrnksTo. Then 
(1231 can be written as 



Aq «sC V47T— <K AdB- 

0-q 



(24) 



The de Broglie wavelength must be long compared to the 
other length scales of the system, in the precise sense defined 
by (|2l. 
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We will adopt the more specific term harmonic scale in- 
variance for this regime, as we will show below that the con- 
dition ( fZ3b restricts the system to a regime where the radial 
trap is well approximated by a quadratic expansion of the po- 
tential about the toroidal minimum. In fact, as we shall show 
in Sec. IIV1 the scale invariant property is more general, and 
the system enters a scale invariant regime whenever V a «: Vq. 



A. Harmonic scale invariance 

A more revealing description of the regime represented by 
d22l i and d23l is obtained by considering the radial curvature 
of the trap: 



<3 2 V(x) 



map;. 



2V ' 



dr 2 

so that at r — r m , we have 
d 2 V(x) 



2r 



2\ 



,- r l<r 



(25) 



> 



dr 2 



— = mwj. 



(26) 







which we have used to define the harmonic radial trapping 
frequency about the minimum of the toroid 



co T = <2^^ = Lo r V21n(y /y (T ). 



(27) 



The toroidal trap potential can be expanded about its radial 
minimum as 

mcozz mcotir - r m ) 
V(x) * V T (x) = V m + + T — (28) 

with error 0{(r - r m ) 3 ). 

Using Eq. d27b for the toroidal frequency, we can write 
Eq. d22b in a more suggestive way: 



r = N (v m -fi)- 



2nr m Ci l2 (S^) 

AdB (PtfitOTWz 



(29) 



where now we see explicitly the scaling with the toroidal 
perimeter 2nr m and the appearance of coj as the physical pa- 
rameter of the radial degrees of freedom. 

Our treatment thus far has relied on the full semiclassical 
expression for the free energy for determining the scale in- 
variant regime. A description of the scale invariant regime 
can also be found by applying the harmonic approximation 
Eq. (l28l l directly to the free energy Eq. (8). In the scale invari- 
ant regime given by Eq. (l23l @{a) can be approximated by the 
quadratic expansion of the potential about the minimum I 
The integral (O, given by 



f 



§{a)= \ dze- am(al * ei ' 1 2n J rire^ w '" 1 

(30) 



can then be approximated as 



0(a) 



In 



ocmcot 



L 



In (y + r m ) dy e 



-amcoj. y~ /2-aV„, 



In 



f 

U —1 



2nr m dye 



-am(j)j y~ j 2—aV m 



2nr„ 



2n 



(31) 



(32) 



(33) 



(34) 



amto z u>T 

Using this with ([8]) gives ( |29l that we previously found us- 
ing the exact semiclassical free energy in the harmonic scale 
invariant regime. 

Having found the scale invariant grand potential through 
two different approaches, we now consider the transition tem- 
perature for the system, which is given by 

_ j_ /^VwjW\ 2/5 / h 2 \ 1/5 

Tc ~k B \ <T(5/2) j \nmrl) 
Defining the toroidal kinetic energy 

na> K = 7. T' 

2nmr~ 1 

and modified geometric mean frequency 

- *i 2 2 ™ 2 2 
(x) = U K (i)j(i)Z = ^0)' r (L),, 

nmo~Q 

the transition temperature becomes 

hcbl N \ 2/5 haj 

T c = — ^0.89 — N 2 ' 5 . (35) 

k B U(5/2)j k B 

This expression closely resembles that for the three dimen- 
sional harmonic trap k B T c = h(co x u) y u> z N / ^(3))' ^ . The scal- 
ing with Af 2 ' 5 is caused by the reduction in the number of 
thermalized degrees of freedom in the system by one, a con- 
sequence of toroidal trapping topology. The characteristic en- 
ergy of the system hCo, which also determines T c , is seen from 
Eq. ( f34b to be scale invariant; invariance is caused by the pre- 
cise dependence of aij on toroidal radius and the fact that this 
degree of freedom is harmonically bound so that it is dou- 
bly weighted in a>. In general harmonically bound degrees of 
freedom have double weight, whereas u>k associated with the 
periodic coordinate is only singly weighted. 



B. Density of states 

The semiclassical density of states is given by 
2t!-(2m) 3/2 



Pie) 



h 3 



JV(£) 



V(x), (36) 



where V*(e) is the spatial volume available to a particle with 
energy e. Making use of cylindrical symmetry, we obtain 



1 m r R+ 
p(e) = T3— drr[e-V eS (r)], 



(37) 
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where V eff (r) = moj 2 r 2 /2 + V <r'' 2/ ^ - V„, ~ mu\{r - r m ) 2 /2 
is the shifted radial potential and energy is now expressed rel- 
ative to V m . The limits of integration are the two solutions of 
£ = V e «(R±), with R < R + . Making the harmonic approxi- 
mation to the potential, the semiclassical density of states be- 
comes 



9(e) 



1 2mr„ 
h 3 u>~ jo 

1 2mr„, 2e 

w 



dy[e — ma>jy 2 /2], 



2e 



(38) 



(39) 



Our analysis of the exact expression of the grand potential 
has shown that this regime is reached rigorously by taking 
Vq <k Va-, and V a <*c k B T simultaneously. Equation ( f39b can 
be written as 



lim 



g(e) = g T (e) = 



4e 3/2 



3 V^(^w) 5 / 2 ' 



(40) 



giving all thermodynamical properties of the system in the 
harmonic scale invariant regime. 



C. Thermodynamics 

Using either (|4D1 . or d34l l and ( [29l , the grand potential can 
now be written as 



T = N (V m -fi) 



ft/2(eW-""') 
PijSHuj) 5 ' 2 



(41) 



For completeness we give the thermodynamic quantities in 
the harmonic scale invariant regime, both above (>) and be- 
low (<) T c . The total atom number in the system is 



N = N + 



(J3hcS) 5 l 2 



(42) 



where iVo > below T c and the condensate fraction is then 
given by 



No =l _ (T\ 512 
N \TJ ' 



again showing the reduction to 5 degrees of freedom. 
The entropy S = -(df/dT)^^ is 



(43) 



Nk B ~ 2 f 5/2 (e^-^)) 



5< 7 £(7/2) / T vV2 
~Nk~B " 2 f(5/2) \l7 



The energy C/ - T + TS + fiN is 

£/> = 5£tM^_^) 
Nk B T 2 {spje^-v,,,)) 



(44) 



(45) 



(46) 



and the heat capacity C = (dU/dT)N,v takes the form 

C> _ 35 ^7/2(^° i ' v '" ) ) 25 &/2^^) 
4 6/2(e^- ,/ '» ) ) 4 f 3/2 (e^- y »)) 



C< _ 35 £(7/2) / T 
Nk B ~ Tf(5/2) 1 7^ 



5/2 



(48) 



(49) 



The discontinuity in the heat capacity across the transition 

AC(T C ) = C>(7/ f ) - C<(7 C ) is 



AC(7V) = 25 £(5/2) 
% 4 £(3/2) 



-3.21. 



(50) 



D. Generalized volume and pressure and equation of state 

A generalization of thermodynamical formalism to trapped 
systems has recently been developed and applied by Romero- 
Rochin et al. l36l l39l l40ll . This formalism provides a def- 
inition of generalized volume and pressure variables for the 
system, allowing an equation of state to be usefully obtained. 
While in principle the formalism allows a generalized volume 
to be defined for any geometry, in general for this system it 
would have to be determined numerically. To gain some in- 
sight into the role of trapping topology in changing effective 
system volume, in this section we apply this approach to treat 
the harmonic scale invariant regime of the ideal gas. 

Proceeding in a similar manner to ll36Tl we see that since 
N, U, S and T are extensive variables and T and fi are inten- 
sive, we may identify the generalized extensive volume in the 
harmonic scale invariant regime as <ir 5 / 2 : 



*y = (u> K u) T (jj z ) 



-1/2 



(51) 



That it takes this form is not so surprising, since at a given 
temperature the gas will be mainly confined to a volume 
of order (toroidal circumference)x(cross sectional area) ~ 
2nr m (k B T/moj 2 z ) 1/2 (k B T/maj 2 T y 12 oc <zr 5/2 . This should be 
compared to the volume for a harmonically trapped gas which 
is of order (k B T /ma> 2 ) 3 ^ 2 oc w~ 3 , where oj = (w,WjW z )'' 3 is 
the usual geometric mean l36ll . We can now identify an impor- 
tant change in the system relative to the purely harmonic case: 
the system has effectively two thermally determined length 
scales parameterizing the cross section, and one purely geo- 
metric length scale that is under direct experimental control, 
the toroidal circumference. 

The conjugate generalized intensive pressure P = 
-(dT I&V)nj is then given by 



V 



(52) 



Above T c this gives the equation of state 'PV = -T as re- 
quired |360. 
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It is important to see that the thermodynamic limit is well 
defined in the scale invariant regime. The appropriate limit is 
Af — > oo, 'V ^ oo, with N/'V = Na> 5/2 — > constant. Since 
'V oc ctqI co r a> z we see that thermodynamic limit requires 



NtL> r io z 



constant, 



(53) 



corresponding to either relaxing the harmonic trapping fre- 
quencies, or increasing the Gaussian beam width cr appro- 
priately. For the latter case, if the system is to remain in the 
scale invariant regime (V lT <K Vo) we must also impose the 
further condition that Vo 

a consistent thermodynamic limit in this regime 



oo at least as fast as cr^ to obtain 



IV. GENERAL SCALE INVARIANCE 

In this section we show that it is possible to further relax 
the condition ( |23l for scale invariance. Our numerical inves- 
tigations show that a scale invariant regime is always reached 
provided V a Vo (see Fig.|2). Thus there are scale invariant 
regimes parameterized by the ratio Va-/k B T ~ 1, correspond- 
ing to A dB ~ a 2 /cr . 

However, the properties of the system are not well de- 
scribed by the harmonic approximation to the potential d28l 
unless k B T «c Vo- also holds. We now seek a general asymp- 
totic expansion of the grand potential which does not require 
this condition to hold. 

We make use of the representation of y(a, x) in terms of the 
confluent hypergeometric function M(a,b,z) l38ll y(a, x) - 
a~ x x"e~ x M{\,\ + a,x), and the asymptotic expansion of 
M(a,b,z) for large real positive z and fixed a,b. To leading 
order in powers of z~ l we easily obtain 

y(a, x) = (ar x Y{a + 1) - f l e*) (1 + 0(z -1 )). (54) 

At leading order in (ySVo) 1 we then find 



k=\ 

Vr 
Vo 



r(kpv a + 1) 

~ (k/3V () ) k l }v - 



(55) 



As Vo/V a -> oo, ( v (^- Vo) ) -> ^- v " } and the second term 
vanishes. We note that, for the vast majority of terms in 
the summation, the inequality 1 «; kfiVa- will hold, even if 
/3Vcr ~ 1. Introducing Stirling's expansion for T{kf5Va- + 1) 
and making use of the identity e^ Vm = e^ v ' r {V cr IV(if v ' J we 
finally obtain the expansion 



lim 



{ v (e^,/3V a ,/3V ) 



^l/2(/ ( "- V '" ) ) 
H l20V<r 

h 288( y 6V (r ) 2 



,(56) 



where the numerical coefficients of the asymptotic expansion 
are the coefficients of the Stirling expansion for T(z). This can 
be used to give a more general expression for the free energy 
than the asymptotic form (l4"TT i which arises from the leading 
term in d56i l. 

In the regime where Vr ^ Vo the grand free energy can 
now be written as 



T = N (V m -fx) 



1 



£>/2(e*"- v ">) 
12fiV<r 



(57) 



This expression for the free energy is our main result and pro- 
vides an asymptotically exact representation of all thermody- 
namics of the system in the scale invariant regime. The energy 
scale Vo only appears through the shift of the potential mini- 
mum V m . Ignoring this trivial shift of energy, we now have 



In^oo ar m 



(58) 



and all ideal gas thermodynamics of the system are scale in- 
variant. In practice this regime is reached rather quickly. In 
this example shown in Fig.|2]the onset of scale invariance oc- 
curs around Vo ~ 3V a . 



A. Non-harmonic corrections to T c 

As an application of the generalized expansion (|57| | we can 
now find a more accurate expression for T c . We easily find the 
following asymptotic expansion for T c : 



1 c 



| k B T c gi/2) (k B T c ) 2 q9/2) 
+ 12V^(5/2) + 288V^(5/2) 



-2/5 



(59) 



where the atom number has been eliminated in terms of the 
harmonic scale invariant transition temperature (f35T > denoted 
by T°. We now expand T c = T° + 8T\ + 6T 2 + which 
we treat formally as a perturbation expansion in powers of 
k B T®/V tT . Solving for the first order correction, we find 



6T[ 



k B r> 



« 1I2 K -1.44X10-^° 



30 V a {(5/2) 



V, r 



(60) 



In Fig.|2]we show the numerical solution of (l20b and compare 
with the harmonic scale invariant form d35l l and the first cor- 
rection for non-harmonic behavior (f60~b . The full numerical 
solution reaches a scale invariant regime for W <k Vo and ap- 
proaches the harmonic behavior for k B T c <K V a . The rapid 
decay of the Stirling expansion coefficients renders the first 
correction given by d60l ) adequate even when k B T° ~ 3^, as 
can be seen from Fig.|2](a) where + 5T}, approaches T c to 
within 3%. 
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FIG. 2: Numerical T c for 87 Rb atoms held in a cylindrically symmet- 
ric harmonic trap for a range of Gaussian potentials parametrized by 
the location of the trap minimum r,„. T c (solid curves) is calculated 
for (a) o"o = 14.1yum, (b) ctq = 24.7/jm, and (c) cro = 49.5//m, by 
solving j20\ . The solid horizontal line in (a) shows T c = 43.3nK 
for the purely harmonic trap (Vb = 0). Dashed lines give the har- 
monic scale invariant T® i35\ , and dash-dotted lines show the cor- 
rection T° + 6Tl given by J601 > . The ratio k B T c (Vo — > <x>)/V a is 
3.5, 0.95, and 0.19 for (a), (b), and (c) respectively, and the verti- 
cal lines show Vo = 3Vo-. Other parameters are N = 10 6 atoms and 
(co-,LO r ) = 2^(15.3,7.8) Hz. 



V. NON-IDEAL CORRECTIONS TO T c 

The ideal gas behavior reveals a scale invariant regime for 
the toroidal trap. While this gives a strong indication of what 
to expect in experiments, it is important to determine how ad- 
ditional effects will change this picture. We now consider the 
two most significant effects on the transition temperature: fi- 
nite size corrections, and interactions. Our main focus here is 
on the modifications these effects will make to the scale in- 
variance of quantum degeneracy. 



A. Mean field interaction shift in T c 



Recent experimental 14111 and theoretical 14211 studies of the 
harmonically trapped Bose gas have established that the value 
for T c is well described by the combination of ideal gas the- 



ory and the first order mean field interaction shift 14311 . Within 
experimental error bars, the mean field interaction shift is the 
dominant effect and critical fluctuations appear to be entirely 
negligible. Since the system we consider is in the semiclassi- 
cal regime we evaluate the first order interaction shift due to 
s-wave collisions from the expression derived by Giorgini et 
al. JH, which is 



8Tf ' 2U J d 3 x dtPjdn [n° (x = 0) - n° (x)] 



1 c 



J dH drPjdT 



(61) 



where «^(x) is the non-interacting semiclassical thermal cloud 
density 



dB k=l 



-Bk( M -V(x)) 
"^72 ' 



(62) 



and Uq = Antra I m gives the interaction strength in terms of 
the s-wave scattering length a. As before, in the harmonic 
scale invariant regime where d23l holds, we can carry out the 
harmonic approximation for the potential. A straightforward 
calculation similar to that of Ref. ll43tl gives 



8Tf aN l/5 ( 8 £(3/2) 2 (l-G) 



a \ 5 V2 7T 1/2 <T(5/2) 6 /5 



(63) 



where a = y/Ti/maj is the modified geometric mean of the 
toroid length scales defined as a 5 = a^a\a\, with aj = 
sjh/mtjj, and 



1 °° 



1 



<T(3/2) 2 f'Wj + k) «3/2) 2 
We can obtain G by writing 

1 



(64) 



j=Uk=l 



-1 t— 1 \ J 



+ k) f'w 2 



aV2Y-S (65) 



to give G - 1/2. Evaluating the numerical factors in 
gives 



ST" 



* r 



-1.53- 



:iN 115 



(66) 



This expression has the same structure as the well known re- 
sult for the shift in T c in the harmonic trap ll43"tl , A different 
characteristic length scale arises here and the dependence is 
now on A^ 5 rather than the A^ 6 behavior in the three dimen- 
sional harmonic trap. As noted in Ref. j43ll . the first order 
interaction shift depends on the geometric mean frequency of 
the trap. Here we see the interaction shift is invariant with the 
size of the toroid, depending on system size only through the 
modified geometric mean length scale a. 



B. Finite size effect on T c 

The effect of finite particle number is calculated to first or- 
der by shifting the ground state chemical potential up to the 



8 



quantum mechanical ground state of the potential 



(67) 



which shifts T c according to 

r f, 



ST. 



1 dT° 



5ji 



com (asm 215 oiT + ah 

5^5/2)1 N J <b 



Evaluating the numerical factors gives 



ST: 



fs 



-0.44 



(68) 



(69) 



As a>T r,„ while a> is invariant, finite size effects can become 
a significant correction for large r m . However, the scaling with 
N~ 2/5 will strongly suppress this effect for large N. 

For the parameters used in Fig.|2](c), the shifts have the val- 
ues: ST' C "'/T° ~ -2 x 1(T 5 (for lOOfim < r m ), and -2 x 1(T 2 < 
ST* S /T° < -6 x 10~ 2 (for lOOjum < r m < 500/jm). In con- 
trast to the harmonically trapped gas, for these parameters the 
finite size correction is the dominant shift in the toroidal trap. 
However, the finite size shift vanishes in the thermodynamic 
limit and for significant atom number (here N = 10 6 ) there 
is a wide range of toroidal radii where scale invariance of T c 
holds to within a few percent. 



VI. REACHING SCALE INVARIANCE 

There are at least three means to reach the scale invariant 
regime with degenerate Bose gases. Firstly, a non-condensed 
gas may be evaporatively cooled below T c into a toroidal 
trap as recently demonstrated experimentally |6J]. A second 
method is to adiabatically load from the harmonic trap into 
the toroid by ramping up the optical potential which we fur- 
ther investigate below. Lastly, a gas in a harmonic trap with a 
small optical potential may be rotated which will push it into 
the scale invariant regime. 

The physics of the rotating case is rather simple: an ideal 
Bose gas in equilibrium in the trap (Q~|i in a frame rotating 
around the z-axis with frequency Q < co r has thermodynamic 
properties of a system in the lab frame with effective radial 
trapping frequency a)± = ^cj 2 . - £2 2 ll44ll . Thus all the proper- 
ties of the system follow from the lab frame analysis presented 
above, after replacing W with an effective rotating frame en- 
ergy V^ = raw 2 cr 2 /2 = Va-(l -Q. 2 /co 2 ). In the toroidal regime 
(where < Vb) the radius to the trap minimum is given 
by 



ln(V /Vcr)+ln 



1/2 



The essential condition for scale invariance becomes 
V„i\ -Cl 2 /co;)« V Q 



(70) 



(71) 



for the rotating gas, and the harmonic approximation will be 
valid when kgT <k V lT (l - Qr/oj 2 ) also holds. The frequency 



1 
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FIG. 3: Adiabatic loading from a harmonic trap into the scale in- 
variant regime, (a) Reduced temperature in the toroidal trap (solid 
line) after adiabatically ramping up the optical potential for a system 
with reduced temperature t H in the harmonic trap. The line t T = t H 
(dashed line) is shown for comparison. The point t H = t* H - 0.67, 
from J76| l. is shown by the vertical line, (b) Condensate fraction in 
the toroidal trap (solid line) after adiabatic loading from initial re- 
duced temperature t H . The condensate fraction for the harmonic trap 
is also shown (dashed line). 



(l34l i determining harmonic scale invariant properties is modi- 
fied to 



n 



to, = 



(72) 



by the rotation. We conclude that, for a gas in rotating equilib- 
rium, scale invariance with respect Vo is preserved, but with 
appropriately modified characteristic energy fta> ± determin- 
ing the density of states. However, since <b ± depends on Q, 
changes in the toroidal radius ( T70b caused by rotation do not 
have the scale invariant property. 



A. Adiabatic loading 

We consider adiabatically loading a degenerate Bose gas 
from the harmonic trap into the toroidal trap. In order to main- 
tain adiabaticity the timescale of loading 7Y should greatly ex- 
ceed the slowest timescale of the system, i.e: min(2^/w ; ) <K 
Tl, Under these conditions entropy will be conserved during 
loading. Equating the entropy in the harmonic trap for the 
regime T < T c 



4<T(4) 3 



s» _ 4<r(4)/r\ J 

Nk B <T(3) \T C ) f(3) lfI ' 
with the scale invariant toroidal entropy of Eq. (l45b 

Si 7f(7/2), 



,5/2 



Nk E 



((5/2) 



(73) 



(74) 
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we find 

/ 8^(5/2) \ 2/5 6/s 
fr(fH) " 17^(3)^(7/2) j (?5) 

for the reduced temperature in the scale invariant toroid after 
adiabatic loading from the harmonic trap. From this expres- 
sion we can find t* H defined by tr{t* H ) = t* H , ie: where adiabatic 
loading has no effect on the reduced temperature in the toroid. 
This is given by 

r* =(7^(7/2)^ 0.67. (76) 

Evaluating the numerical factors gives 

*t = tf /(f H ) 1/S - 1.08&' 5 , (77) 

which is shown in Fig. [3] (a). We see that to a good approx- 
imation tj as tn for the region f# < 1 . Substituting Eq. d77l i 
into Eq. ( 143b . gives the expression 

(^) = l-4/(^) 1/2 -l-l-23^ (78) 

for the condensate fraction in the toroid after loading. This ex- 
pression is shown in Fig.[3](b), where we see that the conden- 
sate fraction is well preserved during loading from the purely 
harmonic trap to the harmonic scale invariant toroidal trap. 

VII. CONCLUSIONS 

We have provided a numerical and analytical treatment 
of the Bose gas trapped in a harmonic-Gaussian trap with 
toroidal topology. We have identified a regime where the 
harmonic-Gaussian toroidal trap has properties that are inde- 
pendent of the radius to the toroidal minimum. In particu- 
lar, the transition temperature to Bose-Einstein condensation 
is scale invariant, and the toroid radius can then be increased 
without altering quantum degeneracy. 

Harmonic scale invariance. — In the regime where k B T c <sc 
Vcr = mco 2 r cr^l2, the radial trap is well approximated by its 
quadratic expansion about the minimum. This regime affords 
an analytical description and we have identified the invariant 



generalized extensive volume *V = (^k^io^Y 111 which de- 
termines the thermodynamics of the system. We have also 
shown that the first order mean field interaction shift to T c 
is scale invariant in this regime. The main limitation on the 
invariance of T c is the finite size shift that vanishes in the ther- 
modynamic limit. 

General scale invariance. — Relaxing the condition k B T c <K 
Vcr, we find that the ideal gas always enters a scale invariant 
regime when V lT <*c Vo, for which the grand potential becomes 
independent of the toroidal radius: dT /dr m = 0. In consider- 
ing corrections to the harmonic approximation to the potential 
we have evaluated the first order perturbation of T c in powers 
of k B T c IVcr, and find it provides an accurate approximation 
even in the regime where V IT ~ k B T c . In practice scale invari- 
ance is reached quite rapidly with increasing Vo, and the onset 
occurrs at about Vo ~ 3V tT for the specific system treated here. 

Reaching scale invariance. — We have considered adiabatic 
loading and applying rotation as ways to reach the scale in- 
variant regime. Adiabatic loading appears to preserve the con- 
densate fraction of the ideal gas quite well, while rotation en- 
hances the height of the Gaussian relative to the rotating frame 
harmonic trap, thus deepening the toroidal potential. Another 
promising method is to evaporatively cool directly into the 
toroidal trap. The existence of a scale invariant regime shows 
that rather than always decreasing with toroidal radius, T c for 
such a system reaches a well defined plateau. A system with 
these properties may be promising for creating and loading a 
large toroidal trap with a persistent current while maintaining 
quantum degeneracy. Future work will focus on more gen- 
eral thermodynamical quantities, the role of interactions, and 
extensions to trapped Fermi gases. 
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